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ABSTRACT 

Stellar winds and other analogous astrophysical flows can be described, to lowest order, by the 
familiar one-dimensional hydrodynamic equations which, being nonlinear, admit in some instances 
discontinuous as well as continuous transonic solutions for identical inner boundary conditions. We 
elaborate on the physics underlying the formation of the discontinuous transonic solutions; discuss 
the characteristics of the time-dependent differential equations of motion to show how a perturbation 
changes profile in time and, under well-defined conditions, develops into a stationary shock discon- 
tinuity; and show how the formation of standing shocks in wind-type astrophysical flows depends on 
the fulfillment of appropriate necessary conditions, which are determined by the conservation of 
mass, momentum and energy across the discontinuity, and certain sufficient conditions, which are 
determined by the flow’s "history". 


L Introduction 


The present paper continues the study of the solution topologies of the time-dependent isother- 
mal wind equations presented by Habbal, Rosner Sc Tsinganos. (1983; Paper I) in these Proceedings 
these two studies are motivated by the desire to understand how the simplest hydrodynamic equatioms 
describi^ spherically symmetric solar-wind type astrophysical flows (Parker 1963) respond to pertur- 
bations which disturi) the radial symmetry of the onitflow (Kopp & Holzer 1976; Munro & Jackson 
1977} or. equivalently (cf. Habbal & Tsinganojs 1983), perturb the outflow by means of a localized 
non-thermal momentum deposition in the thermal gas (Holzer 1977). In this paper, we consider the 
omsequences of a localized phenornienological momentum perturbation term of the form (Paper I) 
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vdiere the free parameters R , o, and t characterize the amplitude, location, spatial and tem- 
pmal width, respectively, of the function D(R,t). The detailed dependence o)f the steady solutions oif 
the governing hydrostatic equations on the parameters D^j, R , and a far T » 0 can be found in 
Habbal &; Tsinganos (1983). The specific purpo>se o)f this noteMS to further elaborate on the physics 
underlying the response of the governing hydrodynamic equations to perturbations of the form (1) 
above (discussed in Paper I), and, in particular, to discuss qualitatively the conditions leading to the 
formation of shocks in the flow. We shall show that in certain limiting cases much of the relatively 
complex behavior can be understood on the basis of the straightforward theory of characteristics. 


2. Nature of the Equations o>f Motion 


b) the following we discuss briefly the nature of the equations of motion and how they allow 
die develojpment o)f shocks in the flow. First, it is convenient to write the continuity and momentum 
balance equations in the following form 
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where c = [dp/dp]^^^ is the isothermal sound speed and A(R,t) and D(R,t) are the given functions of 
the flow tube’s cross section and non-thermal momentum deposition, respectively. 

Eqs. (2) form a system of two quasilinear partial differential equations of the first order for the 
functions p and v of the two independent variables R and t, for which solutions in closed form are 
not known. However, the qualitative behavior of their solutions and the associated physical effects 
can be simply understood by exploiting the hyperbolic nature of these equations; the formal theory 
of hyperbolic partial differential equations shows then that there exist two characteristic directions at 
every point of the (R.t)-plane which define the characteristics of Eqs. (2) (Courant & Friedrichs 1948), 
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The physical significance of these characteristics can be understood as follows. Expressing the 
velocity in terms of the potential $(R,t), 
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we obtain a single differential equation for ^(R,t), 
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Consider then a small perturbation S4(R,t) si^erimposed on some solution $^(R,t) of Eq. (5). 
Following Landau & Lifshitz (1975), we substitute » = + 4$ in (5) and obtain to lowest order the 

following equation for 4$, 

a^4$ a^4$ 2 2 

■ 0 , (6) 

where v^=®^q/ 9R is "background" flow field (see below). We write 4$ = Be , where B is a 
slowly varying function of R and t and the "eikonaT is almost linear in R and t; the group velo- 
city" of the disturbance is then dR/dt = dw/dk, where k = 39ldR and u = d9jdt. We obtain from 

(6) for this "group velocity" of the perturbation, 

dR 

— = Vj,(R,t) ± c , (7) 


a relation which coincides with the slope (3) of the characteristic directions on the plane R-t. From 
the physical point of view Eq. (7) gives the velocity of propagation of sound waves in the moving 
gas relative to some fixed coordinate system; and since a disturbance of finite amplitude and dura- 
tion can be regarded as a superposition of a sequence of small amplitude and duration perturbations, 
Eq. (7) describes the propagation of finite perturbations in the atmosphere as well (Courant & 
Friedrichs 1948, Landau & Lifshitz 1975). It is evident then that the characteristics represent the 
paths of all possible disturbances in the (R-t)-plane. The inevitable result of having a non-uniform 
velocity field, dvJdR * 0, is that the profile of some initial perturbation changes as it propagates in 
the atmosphere. Tor example, if dv^JdR < 0 somewhere in the flow, an outward-propagating wave 
finds its leading segment travelling slower than its trailing edge; the result is that its amplitude pro- 
file "steepens". The perturbation then may "break," Le., it may steepen enough so that the total den- 
sity and velocity are no longer single-valued; in that case the characteristics intersect (Landau & 
Lifshitz 1975), and a shock has been formed in the flow. Whether this indeed leads to an equili- 
brium solution with a discontinuity depends on two facts: first, on the existence of such equilibrium 
solutions involving shocks for the given set of spatial parameters in the steady state equations (Hab- 
bal & Tsinganos 1983); and second, on the temporal history of the flow, as discussed below. 
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3. Evolution of the Perturbation and Shock Formation 


We would now like to see whether the above discussion of characteristics can shed any light on 
an understanding of the behavior of a perturbed isothermal wind flow. We begin with a steady 
and isothermal (single critical point, or Parker) wind flow at t-0, and apply a perturbation of the 
type given in Eq. (1). As shown in Paper I, which of the (maximally three in the present, simplest 
case) distinct steady solutions the perturbed flow finally approaches depends upon the time scale 
peuameter t of the perturbatitm; whereas the character of these possible alternative steady states is 
determined by the remaining parameters of the perturbation [D^ R , and ol which fix the spatial 
form and strength of the perturbation (Habbal Se Tsinganos 198s). How is one of these steady states 
selected by virtue of the time scale of the perturbation? 

The first point to note is that the flow is characterized by three distinct time scales: the advec- 
titm time scale Tj^-R/v, the downstream signal propagation time scale - R/(v+c) [corresponding to 
propagation along the C'*’ characteristics], and the upstream signal propagation time scale t_ ■ R/(v- 
c) [corresponding to propagation along the C~ characteristics]; the crud^ element is that although 
and are comparable in the r^on of interest (i.e, of order of 3 x 10^ sec near the sonic point), t. 
may be very large; the question is then how t OMnpares with these three flow time scales. 


Consider then the evolution of the solutions sketched in Figs. (a)-(d) [see also Fig. 2 of Paper D, 
which qualitatively rqjroduce the essential features of the time-dependent solutions. In each case, we 
have sketched the velocity profiles of the initial Parker-type flow field, v^(R,t-0), the final equili- 
brium flow field V (R,t^) consist«it with some given set of the spatial parameters CD^, Rp, o'] 
characterizing D(R1 smd, finally, an intermediate velodty field v(R,5t) which results after, say, a^time 
it of the order 10’ secs of mom«itum depodtion. The question we would like to address then is: 


given the initial state Vjjj(R,t-0) and the parameters [D^, R , <r, t] of our momentum deposition, 
what is the final state towards which the intermediate state v(R,it) evolves? That is, which of the 


three possible states — the continuous transonic solution, or any one of the allowed (if any) discon- 
tinuous transonic solutions ~ the system chooses to relax to after a suffidently long evolution? 


In the first case (Fig. a), no difficulty arises: in this case, the asymptotic value of D is suffi- 
dently small that the steady perturbed solution (corresponding to t -*”) is unique and continuous, and 
always passes throi^h the outer critical pmnt; the eventual evolution of the intermediate solution is 
hence unambiguous. This is, however, not the case for the three other examples given, for which 
the asymptotic value of D is suffidently large to result asymptotically in three distinct steady tran- 
sonic solutions. In the first of the latter cases (Fig. b), T is large (-40,000 secs) when compared to 
and everywhere, and when compared to T_ almost everywhere. Hence the flow profile evolves 
slowly, and the prqiagation of disturbances along characteristics may be viewed as occurring on a 
background flow field which is relatively steady almost everywhere; the exceptional location is near 
the sonic point(s), where the C characteristic becranes vertical (as can be seen from the sketch of the 
C' characteristics in Fig. f, derived from the numerical solutions given in Paper I). Thus, the flow 
is near equilibrium essentially during the entire course of its history: but which sequence of equili- 
bria is followed? The key point is that this sequence of equilibria (corresponding to a gradually 
increasing amplitude for D^, cf. Habbal & Tsinganos 1983) is characterized by always having an 
outer sonic point downstream from the (eventual) inner critical point, and hence always has a region 
of subsonic flow between these two "x-type" critical points. Hence, once the flow becomes supersonic 
near the (eventual) inner critical point, the C” characteristics are "trapped" between the two critical 
points, and the downstream-fadng C” characteristics emanating just downstream from the inner criti- 
cal point must intersect the upstream-facing C" characteristics emanating from the subsonic region 
just upstream from the outer critical poin^ a shock results (Fig. f), and the subsequent sequence of 
quasi-steady solutions corresponds to the discontinuous solution branch of the steady equations (for 
this choice of T, that corresponding to the upstream shock). Note that from the characteristics one 
may readily obtain the time scale for, viz., fonnation of shocks by locating their first intersection. 


The opposite extreme is obtained for very small T (=10 secs; Fig. d); in this case the time scale 
for the perturbation is far shorter than any signal propagation time scale in the flow. Hence, the 
flow is never near any equilibrium solution during the time of significant change in the perturbing 
momentum deposition term; for 0 < t < T, the flow’s evolution is essentially entirely determined by 
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Figures (a)-(d) are sketches (taken from the solutions obtained in Paper I) of the solution topologies 
for the velocity field of the wind equations (2) with a localized momentum deposition 
function simihr to the g^ussian given by expression (1), for different values of the 
parameters R , 9 and t and at three representative times: t-0, t-6t>0 and r**. 
The dotted smutim is the initial velocity v^(R,t»0), the dashed solution is the inter- 
mediate velocity v(R^) and the solid velocity curves, v (R,r**®), represent the solu- 
tion topologies of the steady state equations (2). All ^ures (aHd) have the same 
values of the parameters R and <r while figures {bHd) correspond to the same set of 
the spatial parameters [D^*R^ er] with T decreasing from (b) to (d). The two verti- 
cal solid lines in figures (oHcR Indicate the location of the steady state shocks. 

Figures (eHf) are sketches of the C characteristics. In Fig. (e) are plotted the fan-like characteristics 
of the initial flow Vj^(R,t-0); and in Fig. (f) the intersecting characteristics of the 
evolving flow field (¥ig. b). The effect of this momentum addition with a large 
value of T is to "turn" downstream the characteristics just downstream of the inner 
aitical point such that they intersect the upstream facing characteristics just upstream 
of the outer critical point 
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the local acceleration due to the imposed perturbation. Thus, by the time the perturbation has 
reached its final amplitude and form, the flow has been locally accelerated (by the external force) to 
virtually its final (steady) spatial form throughout the region of application of the perturbation. 
Note that the discontinuity which temporarily forms is advected downstream from the position of the 
second shock because the C" characteristics immediately downstream from this position face down- 
stream; the continuous transonic solution is obtained. 

Finally, in the intermediate case sketched in Fig. (c), the intermediate solution v(R,5t) finds itself 
to the equilibrium solution which involves the second (downstream) shock. Subsequent propa- 
gation along the C characteristics leads to the steepening of the velocity profile and the intersection 
of the characteristics at a locatim between the two equilibrium shocks; the intersection is convected 
downstream until it finds an equilibrium at the position of the second shock, and the system relaxes 
to the equilibrium solution which involves this second shock. 

In all these cases the discontinuity forms because of the dispersive properties of the non- 
uniform background flow which lead to the gradual intersection of the C characteristics. The value 
of the time constant of the perturbation r determines how soon and where in the flow, in relation to 
the position of the equilibrium shocks, this intersection occurs. In other words, t determines^ how 
rapidly the fan-like C characteristics of the unperturbed flow (Fig. e) will be "turned" so as to face 
one another, Le., converge, and inevitably intersect before the changes in the background flow are 
convected downstream past the positions of the equilibrium shocks (Fig. f). Notice that propagation 
along the C* characteristics does not lead to any interesting effects since the changes are advected 
away relatively fast due to the smaller value of T^. 

4. Suminary 


We have shown that the formation of standing shocks in wind-type astrophysical flows with 
some effective, localized, non-thermal momentum deposition depending on the spatial parameters D^j, 
R and <r - which determine its amplitude, location and width, respectively - and the temporal 
p£umeter t - which determines its temporal width - requires both necessary and sufficient condi- 
tions. The necessary conditions are simply the Rankine-Hugoniot relations for the conservation of 
mass, momentum and- energy - but not entropy - across the shock discontinuity, which in uurn depend 
on the set of the spatial parameters (D*, Rp, The sufficient conditions are determined by the 
detailed temporal evolution of the non-themlal momentum addition and, in particular, on the inter- 
section of the C" characteristics of the How upstream of the position of the equilibrium shocks, 
which in turn depends on the fourth parameter t. It is expected that such standing shocks might be 
formed in astrophysical flows such as the solar and stellar winds, or flows in astrophysical jets, when 
these necessary and sufficient conditions are satisfied. 
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